Andreev Reflection in Ferromagnet/Superconductor/Ferromagnet Double Junction 

Systems 
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We present a theory of Andreev reflection in a ferromagnet/superconductor/ferromagnet double 
junction system. The spin polarized quasiparticles penetrate to the superconductor in the range of 
penetration depth from the interface by the Andreev reflection. When the thickness of the supercon- 
ductor is comparable to or smaller than the penetration depth, the spin polarized quasiparticles pass 
through the superconductor and therefore the electric current depends on the relative orientation 
of magnetizations of the ferromagnets. The dependences of the magnetoresistance on the thickness 
of the superconductor, temperature, the exchange field of the ferromagnets and the height of the 
interfacial barriers are analyzed. Our theory explains recent experimental results well. 

PACS numbers: 72.25.Ba, 75.70.Pa, 74.80.Dm 



I. INTRODUCTION 

The spin-dependent transport through magnetic 
nanostructures has attracted much interest .a In the early 
1970s, Meservey and Tedrow have showed that tunneling 
electrons between a ferromagnetic metal (Fe, Co, Ni) and 
a thin film of superconducting aluminium (Al) are spin- 
polarizedH The ferromagnet /insulator/superconductor 
(FM/I/SC) tunnel junctions are one of the most pow- 
erful tools to extract the spin polarization of the con- 
duction electrons near the Fermi level. In FM/I/SC 
and FM/I/SC/I/FM tunnel junctions, the suppression 
of superconducting gap due to spin accumulation by in- 
jection-^|-st)in- polarized quasiparticles (QP's) has been 
shownBcrmjEllj The QP's spin transport and relaxation 
in SC has been studied in detailBEO 

In recent years, much attention has 
been focused on FM/SG nr - |»GtaUic con- 
tacts botln ^WeticaUyByye-BEl and 
experimentallyllallj'L^LlrEj since the spin polariza- 
tion of conduction-electrons is measured by using the 
Andreev reflectiorc3: An electron injected from the 
FM into the SC is reflected as a hole at the FM/SC 
interface and a Cooper pair is generated in the SC. 
The Andreev reflection includes a conversion process of 
the QP current to the supercurrent carried by Cooper 
pairs in the range of the penetration depth, which is 
approximately eaual to the Ginzburg-Landau, (GL) 
coherence length,EJ from the FM/SC interface.Eil Thus, 
a FM/SC/FM double junction is particularly interesting 
because the magnetoresistance is expected due to the 
overlap of the QP penetration in the SC by the Andreev 
reflection. Recently, Gu et al. have measured the mag- 
netoresistance in a current perpendicular to plane (CPP) 
geometry consisting of a superconducting niobium (Nb) 
thin film sandwiched by ferromagnetic permalloys (Py) 
and proposed a method for estimating ths-penetration 
depth by measuring the magnetoresistanceEj 

In this paper, we present a theory of the Andreev re- 
flection in the FM/SC/FM double junction system and 
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FIG. 1: (a) Schematic diagram of a ferromag- 
net /superconductor/ferromagnet (FM1/SC/FM2) double 
junction system. A superconductor with a thickness of L 
is sandwiched by two semi-infinite ferromagnetic electrodes. 
The system is rectangular and the cross section is a square 
of side W . (b) The current flows along the z-axis. The 
interfaces between FM1/SC and SC/FM2 are located at 
z — —L/2 and z = L/2, respectively. 



derive an expression of the current through the junction 
by extendins-fche theory of Blonder, Tinkham and Klap- 
wijk (BTK).Ea We numerically calculate the current for 
the parallel and anti-parallel alignments of magnetiza- 
tions, and investigate the dependences of the magnetore- 
sistance on the thickness of the SC, temperature, the 
exchange field of FM's and the height of the interfacial 
barriers. It is shown that these dependences are under- 
stood by considering the penetration of quasiparticles to 
the SC by the Andreev reflection process. Finally, we 
compare our.xesults with the recent experimental results 
by Gu et al.H 



II. MODEL AND FORMULATION 

We consider a FM1/SC/FM2 double junction system 
consisting of three rectangular blocks as shown in Figs. 
1(a) and 1(b). The cross section of the system is a square 
of side W and the thickness of the SC is L. The cur- 
rent flows along the z-axis and the interfaces between 
FM1/SC and SC/FM2 are located at z = -L/2 and 



z = L/2, respectively. For simplicity, we assume that 
the system is symmetric: FM1 and FM2 are made of 
the same ferromagnetic materials and the potentials for 
the left and right interfaces are the same. The system 
we consider is described— by the following Bogoliubov-de 
Gennes (BdG) equationcH: 
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where H = — (ft 2 /2m)V 2 — [ip is the single particle 
Hamiltonian, E is the QP energy measured from the 
Fermi energy up and a = + {—) is for the up-(down-)spin 
band. The exchange field h ex (z) is given by 
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where +ho and — ho represent the exchange fields for the 
parallel and anti-parallel alignments, respectively. The 
superconducting gap is expressed as 
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We assume that the temperature dependence 
of the superconducting sap is given by A = 

A tanh (l.74y/T c /T - lj ,El where A is the supercon- 
ducting gap at T — and T c is the superconducting 
critical temperature. In order to capture the essential ef- 
fect of the interfacial scattering, we employ the following 
8- function type potential at the interfaces: 



H{z) = -Z{S{z + L/2) + S{z-L/2)}. (4) 
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Throughout this paper, we neglect the spin-flip scattering 
in the SC and the proximity effect near the interfaces J 2 -!! 

Since the system is rectangular, the wave function in 
the transverse {x and y) directions is given by 



S n i{x,y) ee sin (rmx/W) sin {lwy/W), 



(5) 



where n and I are the quantum numbers which define the 
channel. The eigenvalue of the transverse mode for the 
channel {n,l) is 
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FIG. 2: Schematic diagrams of energy vs. momentum for 
the FM1/SC/FM2 double junction system with the parallel 
and anti-parallel alignments of the magnetizations are shown 
in panels (i) and (ii), respectively. The open circles denote 
holes, the closed circles electrons, and the arrows point in the 
direction of the group velocity. The incident electron with 
up-spin in the channel (n, I) is denoted by 0, along with the 
resulting scattering processes: the Andreev reflection (1), the 
normal reflection (2) at the interface of FM1/SC, the trans- 
mission to the SC (3, 4) and the reflection at the interface of 
SC/FM2 (5, 6), the transmission as an electron to the FM2 
(7) and the one as a hole (8). 



where uq and vq are the coherence factors expressed as 

(8) 
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and k&i ' is the z component of the wave number of an 
electron- (hole-) like QP in the channel (n, I) defined as 

n 2m ~ 
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In the FM region, the solutions are given by 
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where p„ n ? is the z component of the wave number of 
an electron (hole) with a-spin in the channel (n, I); 

V2777 



Pa,nl 



'■ a/mf ± E ± ah ex - E n 



(11) 



The solution of the BdG equation (|l|) in the SC region 
is given by 
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The wave function of the FM1/SC/FM2 double junction 
system is given by the linear combination of the solutions. 
Let us consider the scattering of an electron with up-spin 
in the channel {n, I) injected into the SC from the FM1 
(0 in Figg). There are the following eight processes: the 
Andreev reflection (1 in Fig.||), the normal reflection (2 
in Fig.0) at the interface of FM1/SC, the transmission to 



3 



the SC (3, 4 in Fig.^) and the reflection at the interface of 
SC/FM2 (5, 6 in FigM), the transmission as an electron 
to the FM2 (7 in Fig.0) and the one as a hole (8 in Fig.|). 
Therefore, the wave function in the FM1 (z < —L/2) is 
given by 



0; 



S n i(x,y). (12) 



In the SC {-L/2 < z < L/2) we have 



and in the FM2 (L/2 < z) 



Sulfa v), 



(13) 
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Here p^„ ; , ^ and g* nZ are the wave numbers in 
the FM1, SC and FM2,' respectively. The coefficients 
Q>(j,ni7 ^o",ni) C(7,ni 7 ^(r,nl 7 &<r,nl •> Po",ni 7 ^cr.nh ^md T\ a n \ are de- 
termined by matching the wave functions at the left and 
right interfaces. The matching conditions for the wave 
functions ( |l2| ) - (p|) at the interfaces are 
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Solving Eq.(|l5|), the probabilities of transmission-and re- 
flection are calculated following the BTK theoryEj When 
an electron with a-spin is injected from the FM1, the 
probability of the Andreev reflection R^ e n i, the normal 
reflection R ee , and the transmission as an electron and 

<r,ru 

the one as a hole to the FM2, T*^ and are given 



by 



Pa.nl 
Pa.nl 



R 



(E) 



= b 



rpe e 
ff.nl 



Ta'nl (E) 



Pcr,nl 



dt 



il c a,nli 



(16) 



where the subscript e'(h') in Eq. ( |l6| ) indicates the elec- 
tron (hole) in the FM2. Let us evaluate the current in 
the FM1. When the bias voltage V is applied to the 
FM1/SC/FM2 system, the current carried by electrons 
with (7-spin is written as 

/>oo 

£ = [U(E)-f^(E)]dE, (17) 

nl J ° 

where h is Planck constant and /_» (E) is the distribution 
function of an electron with a positive group velocity in 
the z direction and expressed as 

/-. (E) =fo(E-s£), (18) 

where fo (E) is the Fermi distribution function. The dis- 
tribution function of the electron with a negative group 
velocity in the z direction is written as 
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where v a nl is the Fermi velocity of an electron with a- 
spin in the channel (n, I) of the FM1 (FM2) and T>^ is 
the density of states of er-spin band in the channel (n, I) 
of the FM1 (FM2). Using the conservation of probability, 
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The current carried by holes is calculated in the similar 
way. 

The total current in the FM1/SC /FM2 double junction 
system is obtained as 
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FIG. 3: MR as a function of the thickness of the SC, kpL. 
From top to bottom, temperature T/T c is 1, 0.9, 0.7, 0.5, 0.3, 
and 0.1. We assume £q(E = T = 0) = 200/k F . 



Note that this expression of the current, Eq. (^l|) re- 
duces to the one derived by Lambertca for the nor- 
mal metal/superconductor/normal metal system when 
h = 0. 

The magnetoresistance (MR) is defined as 



MR: 



Rap — Rp 
Rp ' 



(22) 



where i?p(AP) = V/Ip(AP) is the resistance in the parallel 
(anti-parallel) alignment. 



III. RESULTS 

In Fig. H the MR is plotted as a function of the thick- 
ness of the SC, L, multiplied by the Fermi wave number 
kp. We assume that the strength of the interfacial bar- 
rier Z = and the exchange field ho = 0.5/xf- The side 
length of the cross section is taken to be W = 1000//cf 
throughout this paper. When the SC is in the normal 
conducting state (T/T c = 1), the MR is constant since 
we neglect the spin-flip scattering in the SC. When the 
SC is in the superconducting state (T/T c = 0.1, 0.3, 0.5, 
0.7, and 0.9), the MR decreases with increasing the thick- 
ness of the SC. The MR at low temperatures T/T c < 1 
shows an exponential decrease in a wide range of fcpL. 
The decrease of the MR due to the superconductivity can 
be explained by considering the decay of the quasiparticle 
current in the SC. To obtain the charge and spin currpits, 
we extend the BTK theory of the Andreev reflections to 
the FM1/SC/FM2 double junction system. 

Let us first consider the charge transport in the SC. 
When the current flows in the positive z direction, elec- 
trons and holes are injected into the SC from the FM1 
and the FM2, respectively. The conservation law of the 
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FIG. 4: Spatial variation of the z component of (a) the charge 
current density and (b) the spin current density in the SC 
with the thickness L = 3000/fcf is shown. jc^ AP \ ipi^' 

and j^* P ' are the QP current density, the supercurrent den- 
sity and the spin current density in the parallel (anti-parallel) 
alignment, respectively. We assume £q(E = T = 0) = 
200/fc F . 



charge density Q a = e£ ni (|/£| 2 - |«£ ? | 2 ) in the SC, 
where f° d and are electron- and hole-like components 
of the wave function in the channel (n, l,cr), respectively, 
is derived from the BdG equation (|l|) and obtained as 

dQ ° - " yEM/~V^). (23) 
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where Jq is the QP current density with cr-spin and the 
z component of Jq per unit area jq is written as 
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The right hand side of Eq. ( |23|) corresponds to the gra- 
dient of supercurrent carried by Cooper pair Jp air defined 

as 
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FIG. 5: The temperature dependence of the QP current pen- 
etration depth £q for E — 0. We assume (,q{E = T — 0) = 
200/& F . 
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The z coordinate dependences of the QP current den- 
sity JQ = 3q+3q and the supercurrent j pair = J^r+ipair 
in the case that the thickness of the SC L = 3000/fci? are 
shown in Fig. |](a). We find that j'q decays from the 
interfaces of the FM1/SC and the SC/FM2 and becomes 
zero in the interior of the SC. On the other hand, j pa ir 
increases and becomes dominant in the interior of the SC 
to conserve the total current density. In the energy re- 
gion below the superconducting gap (E < A) where the 
energy of the transverse mode E n i is smaller than Fermi 
energy fxp, the wave number k nl is expanded as 
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The imaginary part in Eq. ( p7| ) gives the exponential de- 
cay term exp(— z/£q) in j'q, where £q is the penetration 
depth given by 



2VA 2 - E 2 ' 



(28) 



where vp is the Fermi velocity. £q has a strong tem- 
perature dependence shown in Fig. || When the thick- 
ness of the SC is much larger than the penetration depth 
(L » £q) as in Fig. 0(a), Jq decays in the range of £q 



from the interfaces. Note that £q is approximately equal 
to the clean-limit GL coherence length £(T) in the low 
energy regime: £ Q (E = 0) ~ 1.2£(T)£3 

Next, we consider the spin transport in the SC. The 
conservation law of the spin density S = Pf — Py , where 

P<j =Eni(\fni\ 2 + l5n/| 2 )- & derived from the BdG equa- 
tion (jl]) and expressed as 
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where J sp m = Jp — Jp is the spin current density. The 
z component of J P per unit area jp is written as 
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Figure f|(b) shows the z coordinate dependence of the 

T I 

spin current density j S pin = jp — jp in the SC with the 
thickness L = 3000/fcF. As shown in Fig. |](b), no spin 
current flows through the SC both in the parallel and in 
the anti-parallel alignments because the QP current with 
spin changes to the supercurrent carried by Cooper pairs 
with no spin. This means that the spin injected from 
FM1 does not reach to the FM2. As a result, the QP 
current density in the parallel alignment jg and that in 
the anti-parallel alignment jg P are almost the same as 
shown in Fig. |^(a). 

The z coordinate dependence of the charge current in 
the case that the thickness of the SC is much smaller than 
the penetration depth (L <g; £q) is shown in Fig. ||(a), 
where the QP current density j'q is almost constant and 
the supercurrent density j pa i r is nearly zero in the SC. As 
shown in Fig. ^(b), the value of the spin current in the 
parallel alignment is larger than that in the anti-parallel 
alignment because the value of QP current with up-spin 
is much larger than that with down-spin in the parallel 
alignment, whereas the value of QP current with up-spin 
is equal to that with down-spin in the anti-parallel align- 
ment. This means that the spin injected from the FM1 
is transferred to the FM2 and therefore the value of the 
QP current density strongly depends on the relative ori- 
entation of the FM's magnetizations shown in Fig. ||(a). 

From the above discussion, the result shown in Fig. 
H is understood as follows. In the SC, the QP current 
with spin decreases exponentially and changes to the su- 
percurrent carried by Cooper pairs with no spin in the 
range of £q from the interfaces. As a result, it becomes 
difficult to transfer the spin from the FM1 to the FM2 
and the MR decreases with increasing the thickness of 
the SC. The finite MR in the region of large L is due 
to the QP's with energy above the superconducting gap 
(E > A). 

The temperature dependences of the MR normalized 
by the value at T c (MR norm ) for several values of the 
exchange field ho in the case of L = 3000/fcp and 
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FIG. 6: Spatial variation of the 2 component of (a) the charge 
current density and (b) the spin current density in the SC with 
the thickness L = 10/&F is shown. The parameters are the 
Fig. g 



same as m 



L = 300/kp are shown in Figs. 0(a) and 0(b), respec- 
tively. The MRnorm decreases with decreasing tempera- 
ture because the number of electrons and holes with the 
energy E > A which contribute to the MR norm decreases 
with decreasing temperature. At low temperatures, elec- 
trons and holes mainly distribute in the energy region 
E < A. When the thickness of the SC is much larger 
than the penetration depth (Fig. 0(a)), electrons and 
holes with energy E < A injected to the SC do not con- 
tribute to the MRnorm because the QP current changes 
to the supercurrent in the SC by the Andreev reflection 
and the QP transmission from the FM1 to the FM2 does 
not occur. As a result, the MR norm becomes zero at 
low temperatures T/T c < 0.4. On the other hand, when 
the thickness of the SC is comparable to the penetration 
depth (Fig. 0(b)), the QP transmission by the Andreev 
reflection in the energy region E < A occurs and there- 
fore the finite MR norm remains even at low temperatures. 

Figures ^(a)-(c) show the temperature dependence of 
the MR in the cases of Z = 0, 1 and 3, respectively, for 
the several values of L. The temperature dependence of 
the MR in the case of the transparent interfacial barrier 
(Fig. |^(a)) is explained by the same way as in Fig. 0. In 



FIG. 7: MR normalized by the value at T c (MR norm ) is plotted 
as a function of temperature. The thickness of the SC is taken 
to be IzfL — 3000 and 300 in panels (a) and (b), respectively. 
From bottom to top, the exchange field ho/fiF = 0.1, 0.3, 0.5, 
0.7, and 0.9. The MR at T = T c is plotted as a function of 
the exchange field ho in the inset of panel (a). We assume 
£ Q (E = T = 0) = 200/fc F . 



the case of the finite interfacial barrier (Figs. ||(b) and 
^(c)), for all values of L, the normal reflection mainly 
occurs especially in the energy region below the super- 
conducting gap (E < A) because of the scattering at the 
interfaces. Therefore, the main contribution to the MR 
at temperature 0.3 < T/T c < 1 comes from the QP's 
transmission in the energy region E > A, whose proba- 
bility is independent of L. As a result, the differences in 
the magnitude of the MR for the different values of L be- 
come smaller especially for temperature 0.3 < T/T c < 1. 



IV. COMPARISON WITH EXPERIMENT 

Let us compare our theory with recent experimental 
results in Py/Nb/Py structures measured by GLu et al..c3 
The mean free path in the Nb film I ~ 6 nmEj is much 
smaller than the clean-limit coherence length £0 ~ 40 
nniE3 and therefore the Nb film is in the diffusive regime. 
In order to analyze the experimental results in the dirty 
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FIG. 8: MR for the system with a several value of interfacial 
barrier is plotted as a function of temperature T/T c . The 
strength of the interfacial barrier is Z = 0, 1 and 3 for panels 
(a), (b) and (c) respectively. From top to bottom the thick- 
ness of the SC is k F L = 300, 600, 900, 1200. We assume 
£ Q {E = T = 0) = 200/fc F . 



Nb film, we need to extend the theory in the ballistic 
case to that in the diffusive case. The diffusive effect on 
the Andreev reflection is incorporated into our theory by 
replacing the penetration depth £q in the ballistic theory 
with the penetration depth in the dirty-limit £q. Thus, 
the value of £q(E = 0) at T = obtained by fitting 




T/T 



FIG. 9: MR is plotted as a function of temperature T/T c . The 
solid curves show theoretical results for the thickness of the 
SC, L =30, 40, 50, 60, 80, and 100 nm from top to bottom, 
where Uf is taken to be 1 A -1 for iNb. The symbols show the 
experimental results by Gu et ai£3 for the thickness of the 
Nb, t Nb =30, 40, 50, 60, 80, and 100 nm from top to bottom. 



the experimental data is interpreted as the dirty-limit 
penetration depth £g(E = 0) - l.2yJTfa£ Q (E = 0)Eil 
at T = 0. Figure pTshows the excess resistance AR = 
Rap — Rp normalized by the value at (in the experiment, 
T slightly above) T c (Ai? norm ) as a function of tempera- 
ture. The solid curves indicate the calculated-results and 
the symbols indicate the experimental ones £3 By fitting 
the calculated values to those of the experimental data, 
we obtain £%(E = T = 0) = 46, 36, 36, 33, 30, and 27 
nm for the curves of L=30, 40, 50, 60, 80, and 100 nm. 
respectively, where kp is taken to be lA _1 for Nb.El 
The values of the penetration depth £q(E = T = 0) 
estimated by our theory become larger for the smaller 
thickness of the Nb film and are larger than the dirty- 
limit penetration depth in a bulk Nb ~ \/£xil = 16.2 nm. 
This indicates that A in the Nb film is reduced com- 
pared to that in a bulk Nb. The suppression of A is due 
to the proximity effect. Actually, the height of the real- 
istic superconducting gap depends on the position z in 
the Nb film by the proximity effect. Here we interpret 
the value of A as the averaged value of the realistic su- 
perconducting gap with respect to z in the Nb film. Gu 
et al. have obtained the dirty-limit penetration depth 
£q{E = T = 0) ~ 16.5 nm by assuming that Ai? norm 

is proportional to exp(— tNb/£,Q ), where t^b is the thick- 
ness of the Nb. This value of the penetration depth is 
comparable to those estimated by our theory. 

Although we neglect the effect of spin relaxation on 
the MR in our theory, this assumption is justified as fol- 
lows. In the Nb film in the Py/Nb/E* structure, the 
spin diffusion length X s is about 50 nm,c3 while the pen- 
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etration depth £q(T = 0) is about 27 ~ 46 nm at low 
energy. As seen in Fig. ||, £q(T) is almost constant at 
low temperatures and shows the divergent behavior only 
near T c , indicating that £q is smaller than A s in most of 
the temperature range below T c , Therefore, the effect of 
the QP current penetration is dominated for the MR and 
the spin relaxation effect on the MR is neglected except 
in the close vicinity to T = T c . 

V. CONCLUSION 

The magnetoresistance in the ferromag- 
net/superconductor/ferromagnet double junction 
system is studied theoretically. The dependences of the 



magnetoresistance on the thickness of the superconduc- 
tor, temperature, the exchange field of ferromagnets 
and the height of the interfacial barriers are understood 
by considering the Andreev reflection of spin-polarized 
current. Our theory shows good agreement with recent 
experimental results. 
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